Immersed boundary method for high-order flux reconstruction based on volume penalization

In the last decade, there has been a lot of interest in developing high-order methods as viable option for unsteady scale-resolving-simulations which are increasingly important industrial design process. High-order offer advantages low numerical dissipation, high efficiency on modern architectures and quasi mesh-independence. Despite significant advance solution methods, general CFD workflow (geometry, CAD preparation, meshing, solution, post-processing) largely remained unchanged, with mesh generation being bottleneck often determining overall quality solution. this work, we aim to combine Flux Reconstruction (FR) method simplicity (or lack thereof) Immersed Boundary Method (IBM) steady problems over moving geometries. The volume penalization (penalty-IBM) is selected its ease implementation robustness. Detailed discussions about implementation, including boundary representation, mask function, data reconstruction, selection parameter given. Advantages combining framework shown by various test cases. approach firstly demonstrated linear advection-diffusion equation investigating convergence coupled FR-IBM approach. Thereafter, accuracy canonical (static) cases 2D 3D when compared standard body-fitted unstructured simulation. Finally, handle geometries flow around an airfoil pitching plunging motions.